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ABSTRACT 

We describe a new method for robustly testing theoretical predictions of red giant evolution near 
the tip of the giant branch. When theoretical cumulative luminosity functions are shifted to align 
the tip in /-band and normalized at a luminosity level slightly brighter than the red giant bump, 
virtually all dependence on age and composition (heavy elements and helium abundance) is eliminated. 
While significant comparisons with observations require large samples of giant stars, such samples are 
available for some of the most massive Milky Way globular clusters. We present comparisons with 
the clusters NGC 2808 and M5, and find that NGC 2808 has a deficiency of bright giants (with a 
probability of less than about 3% that a more extreme distribution of giant stars would have happened 
by chance) . We discuss the possibilities that underestimated neutrino losses or strong mass loss could 
be responsible for the deficit of giants. While we cannot rule out the neutrino hypothesis, it cannot 
explain the apparent agreement between the M5 observations and models. On the other hand, strong 
mass loss provides a potential link between the giant star observations and NGC 2808 's unusually 
blue horizontal branch. If the mass loss hypothesis is true, there is likely a significant population of 
He white dwarfs that could be uncovered with slightly deeper UV observations of the cluster. 
Subject headings: neutrinos — stars: evolution — stars: luminosity function — stars: mass loss — 
globular clusters: individual (NGC 2808, M5) 



1. INTRODUCTION 

At one level, our understanding of the late stages of 
the evolution of a low-mass giant are solid: increasing 
electron degeneracy in the core coupled with increasing 
temperature results in the "fiash" ignition of helium, ter- 
minating the red giant branch (RGB). However, some 
of the physics inp uts to the models retain some signifi- 
cant uncertainties (jBiork fc Chaboveil2006l : ISalaris etall 
I2OO2I ). These uncertainties (such as electron conduction, 
neutrino emission, and mass loss) have significant effects 
on the observable characteristics of the brightest stars in 
old stellar populations, with corresponding influences on 
population synthesis models. Testing stellar models for 
the late RGB is difficult though because the stars evolve 
on short timescales make up a small fraction of all of the 
stars in a population. 

Because NGC 28 08 is a very massive globular clus- 
ter (My = —9.39; iHarrisI [T996f ) and because the stars 
within it appear to have nearly uniform metal content, 
it provides us with one of the largest "clean" samples 
of stars for examining the late phases in stellar evolu- 
tion. At the same time, NGC 2808 has an extremely 
peculiar bimodal distribution of horizontal branch (HB) 
stars that has been known since the fi rst moderately- 
deep photometry was taken (|Harrislll974[ ). The first deep 
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color-magnitude diagram (jSosin et al.lll997| ) was an even 
greater shock, revealing a blue HB tail extending 3.5 mag 
fainter in V and containing two additional gaps in the dis- 
tribution of stars. As such, NGC 2808 is one of the more 
obnoxious examples of the "second parameter" problem, 
in which HB star distributions cannot be explained based 
on metallicity alone. While some aspects of the HB dis- 
tribution have been explained since the original observa- 
tions, the overall bimodality and extended blue HB tail 
have not. The most recent ph otometric study of the clus- 
ter bv lCastellani et aP (|2006f ) used archival Hubble Space 
Telescope (HST) images along with wide-field ground- 
based observations. Their sample can be improved upon, 
however, because they only used a fraction of the HST 
observations available. 

2. OBSERVATIONS AND DATA REDUCTION 

Three HST datasets were used in this work. The only 
dataset that had not been previously published was ob- 
tained with the HRC detector of the Advanced Camera 
for Surveys (ACS) using FiSbW and F555W filters as 
part of proposal ID 10335 (P.I. H. Ford). We also re- 
duced archival WFPC2 images from proposals 6095 and 
6804 (Pis Djorgovski and Fusi Pecci). The HRC images 
mostly overlap the fields of the PC chips, but allowed 
us to resolve a number of stars that were blended in the 
WFPC2 images. 

The HRC frames were processed using the DOLPHOT 
photometry package ^ with its module tuned for ACS 
data. Individual frames (prior to drizzling in the ACS 

^ http:/ /purcell. as. arizona.edu/~andy/dolphot/ 
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pipeline) were obtained from the ACS team in order to 
get photometry on the giant stars. For the photometry, 
differences in effective pixel area were corrected within 
DOLPHOT through the use of a pixel area map pro- 
vided on the ACS website ^. For relative astrometry, 
we used the fourth-order polynomial corrections pro- 
vided in the most recent IDCTAB file for the dataset 
(q692007bjJdc.fits) to facilitate the creation of a com- 
mon coordinate system with the WFPC2 images. The 
WFPC2 images were analyzed using the HSTPhot pho- 
tometry package Images from proposal 6095 were 
taken in the F218W, FA39W, and i^555W^ filters, while 
those from the proposal 6804 were taken in the F160BW, 
i^336W^, i^555W^, and F814W^ filters. F439W^ and 
F336W observations were particularly useful for separat- 
ing AGB stars from the RGB. Known geometric distor- 
tions in WFPC2 images were corrected using the MET- 
RIC task in the STSDAS package within IRAF. 

Of the filters available, i^814W^ is most appropriate 
for giant branch study since near-infrared filters are 
most directly correlated with luminosity. To maximize 
the size of our sample we used common stars from the 
two WFPC2 fields to derive a transformation from the 
(?7i439— JTisss, TO555) color-magnitude diagram to msu us- 
ing a second-order polynomial containing magnitude and 
color terms. The residuals show no systematic trends 
with magnitude or color over the range considered. The 
median residual was —0.001 mag for the 814 stars in com- 
mon, with a semi-interquartile range (half of the differ- 
ence between the 25th and 75th percentile in the ordered 
list of residuals) of 0.030 mag. 

3. ANALYSIS 

The very tip in / has been used as a distance indi- 
cator for large, old stellar populations largely thanks to 
its very weak dependence on age and ch emical composi- 
tion (see lDa Costa fc Armandrofilll990l for an early ref- 
erence). This can be seen empirically in CMDs of com- 
posite stellar p opulations i n nea r by galaxies (for recent 
examples, see iRizzi et al.l 120061 : iMouhcine et al.l 120051 : 
iBellazzini et al.ll2005l l. Theoretical predictions confirm 
that t he composition depende nce is very weak for [Fe/H] 
< -1 (|VandenBerg et al.l[2006( ). 

Even with the large number of RGB stars in NGC 2808, 
it is worthwhile to find the most robust comparison with 
models possible. The tests below involve cumulative lu- 
minosity functions (LFs), counting stars starting at the 
observed TRGB. A comparison with models requires a 
magnitude shift (akin to the distance modulus) and a 
vertical normalization. For the horizontal (magnitude) 
shift, we have opted to shift the models to the observed 
magnitude of the brightest cluster giant. As noted above, 
this has the advantage that it is largely independent of 
chemical composition and age (although its absolute po- 
sition is affected by uncertainties in physics that affect 
the timing of the helium flash). For the vertical normal- 
ization, we have forced the models to have a number of 
stars equal to the number of observed stars just brighter 
than the RGB bump at msi4 ~ 14.8. As seen in Fig. [U 
by normalizing in this way, the models for a wide range 
of compositions (heavy elements and helium) and ages 

^ http:/ /www. stsci.edu/hst/acs/analysis/PAMS 
^ http:/ /purccU. as. arizona.edu/~andy/hstphot/ 



overlie each other to a high degree. To state this an- 
other way, the late evolution of RGB stars is virtually 
independent of input parameters, and results from the 
strong correlation between luminosity and mass of the 
helium core. 

While there is little dependence on age or composi- 
tion, there are differences in physics from model set to 
model set that affect the cumulative LF. In the bot- 
tom right panel of Fig. [H we com pare the predic- 
tions from models by the Yale-Yo nsei ([Demaraue et al.l 
120041) ■ Teramo (.Cassisi et al.l l2004D. and Victoria-Regina 
(jVandenBerg et al .1 12006^ models. The Teramo and 
Victoria-Regina predictions overlie each other almost 
perfectly, while the Yale-Yonsei models predict consid- 
erably more giants (and therefore slower evolution) near 
the TRGB. The two Yale-Yonsei models use different 
bolometric corrections. While this does affect the cu- 
mulative LF, it does not account for the entire differ- 
ence with the Victoria-Regina and Teramo models. The 
implementation of neutrino energy loss rates is proba- 
bly respons ible for the differe nce: the Yale-Yonsei mod- 
els employ lltoh et aTl ( 19891) rates, while the Teramo 



models use (jHaft et al.l 119941) rates and the Victoria- 
Regina models use lltoh et all (|19 96 ) rates. The dom- 
inant plasma neutrino emissi on rates have b een up- 
dated several times since the lltoh et al.l (|1989t ) paper, 
and the values used in the Yale-Yonsei isochrones are 
probably too low, which allows giants to evolve more 
slowly because nuclear reactions don't need to provide 
for larger neutrino energy losses. As corroboration, we 
note that the Yale-Yonsei models have lower luminosi- 
ties at the TRGB (log(L/L0) « 3.33) than the Teramo 
and Victoria-Regina models (log(L/L0) 3.37), consis- 
tent with lower cooling rates and an earlier flash. This 
is an important issue because systematic uncertainties in 
the characteristics of the helium flash translate to cor- 
responding uncertainties in the properties of HB stars, 
including their luminosities. 

In Fig. [21 we compare the observatio ns for NGC 2808 
and for M5 (jSandquist fc Boltd [200I to the Victoria- 
Regina models (though the results are nearly the same 
for the Teramo models). To judge the significance of the 
differences, we used two methods to estimate the prob- 
ability of detecting a smaller sample of bright giants if 
they were drawn from a cumulative distribution given by 
the Victoria-Regina models — these are essentially esti- 
mates of the probability of a "false alarm". First, we 
did one-sided Kolmogorov-Smirnov (K-S) tests at vari- 
ous brightness levels. This was necessary because ad- 
justment of the faint limit for the giant sample affects 
the size of an absolute deviation in the cumulative dis- 
tributions (the statistical quantity used in the test) as 
well as the total sample size, both of which affect the 
computed significance of a deviation. When the faint 
limit is set 1.8 magnitudes below the TRGB, the RGB 
sample is 142 stars for NGC 2808, the maximum abso- 
lute deviation is 0.11, and the probability that a sample 
of the same size with a more extreme deviation would be 
drawn from the theoretical distribution is 4.3%. By com- 
parison, for M5 the sample was 91 stars, the maximum 
absolute deviation was 0.06, and the probability was 89% 
that a sample with a more extreme deviation would be 
drawn from the theoretical distribution. By comparison, 
there is a 0.9% chance that a distribution drawn from the 
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Yale-Yonsei model would deviate to a larger degree than 
the observed NGC 2808 distribution, and a 28% chance 
that one would deviate more than the M5 distribution. 

Our second method used predictions from a binomial 
distribution function. Salaris & Cassisi' (1997) used a 
similar formulation to determine the likelihood of finding 
stars within a certain magnitude range near the TRGB. 
If we have a total sample of RGB stars Nrcb brighter 
than a certain level, then the theoretical cumulative lu- 
minosity function can be used to predict the probability 
Pi that any one star in that sample will be found in the 
brighter portion of the sample (brighter than a second 
faint limit that is closer to the TRGB). If the theoreti- 
cal cumulative luminosity function provides an accurate 
model of the relative evolutionary timescales for the gi- 
ants, then the probability of measuring a number of stars 
n < in the brighter portion of the sample is 



P. 



<N 



N 

E 



RGB'- 



^ n\{NRGB - n)\ 



:Pni-p, 



For NGC 2808, we used the RGB sample between the 
TRGB and mgu = 14.8 just brighter than the RGB 
bump {Nrgb = 441). We find that the probability 
reaches a minimum of about 1.3% for a faint limit of 
771814 = 12.4 {N = 36 and Pi = 0.115), but is less than 
about 5% for msu < 12.5 {N = 46). Both of our meth- 
ods agree that the probability of a false alarm for NGC 
2808 is (conservatively) less than 5%. 

The horizontal (magnitude) shift used to align the 
model and observed TRGB has a small effect on the 
calculated probabilities. Because we use the brightest 
observed giant, we underestimate the luminosity of the 
TRGB by some amount. However, because the models 
appear to be overestimating the number of bright giants, 
this only ends up making our test probabilities conser- 
vative overestimates. The probability of finding at least 
one star within a certain magnitude of the TRGB can 
be determined from the binomial distri bution and an as- 
sumed cumulative luminosity function (jSalaris fc Cassisil 
[1997). Using the Victoria- Regina models, we find there 
would be a 50% probability of having at least one star 
within 0.02 mag of the TRGB for a sample of the same 
size as the one in NGC 2808. If the brightest NGC 2808 
giant is as much as 0.05 mag fainter than the TRGB (the 
binomial distribution predicts less than a 5% chance of 
this), the K-S probability would be reduced to 2.6%. The 
probabilities derived from the binomial distribution in- 
crease slightly if the TRGB brightness is underestimated: 
a 0.05 mag underestimate increases the minimum prob- 
ability from 1.3% to 2.1%. Thus, the deficit of bright 
RGB stars in NGC 2808 seems secure. 

We have two possible explanations for the observa- 
tions. First, larger neutrino emission rates would ac- 
celerate evolution near the TRGB, as mentioned above 
with regard to the Yale-Yonsei models. The NGC 2808 
observations may be showing this clearly, while statis- 
tical fluctuations may be concealing it in M5. We have 
no reason to believe that RGB neutrino emissions should 
differ systematically from cluster to cluster, so this ex- 
planation would disconnect the deficit of RGB stars from 
discussions of HB morphology. The /-band location of 
the discrepancy in the NGC 2808 CLF (starting 1.6 
mag below the TRGB) holds some physical information 



if it results from greater-than-predicted cooling of giant 
cores. As long as the temperature and density sensitiv- 
ities of the plasma neutr ino emission rates agree with 
those of lHaft et al.l ()1994f ). modifications of the emission 
rate are unlikely to explain the NGC 2808 observations 
because they would primarily affect RGB stars within 
about 1 mag of the TRGB. For typical RGB core con- 
ditions, the energy loss rate has fairly large temperature 
and density depend ences (g oc p^-^T^, evaluated from 
the lHaft et al.|[l994l formulae) that are responsible for a 
rapid increase in energy loss. Of course, if there is an 
energy loss mechanism having different density and tem- 
perature dependences, then the shape of the LF would 
be changed. While we cannot rule out the possibility of 
systematic errors in the theoretical neutrino loss rates 
used in the stellar models, new physics would probably 
be required to match the observations. 

A second explanation involves enhanced mass loss near 
the TRGB for a fraction of the stars. Models predict that 
stars can leave the RGB before He flash if the mass of the 
star's envelope decreases below a critical value. These 
stars can have a "hot flash" He ignition after leaving the 
RGB, before finally appearing on the blue end of the HB 
for a time approximately as long as that of a typical HB 
star. With even more mass loss on the RGB, a star will 
leave the RGB earlier, a hot He flash will be prevented, 
and a He white dwa rf would result. According to models 
(|D'Cruz et al.lll996D . stars leaving the RGB before they 
get within 0.4 mag of the TRGB will produce He white 
dwarfs. In NGC 2808, the vast majority of the stars seen 
on the HB in NGC 2808 should have gone through the 
entire RGB phase (complete with He flash) because their 
positions on the HB imply substantial envelope mass. 
However, the deflcit of stars more than 0.4 mag fainter 
than the TRGB would imply we should expect to have a 
significant population of He core white dwarfs if the mass 
loss hypothesis is true. Based on the lack of a blue HB 
tail in M5, we would not expect an RGB deficit there. 

The population ratio R = Nhb/Nrgb = 1.62 ± 0.07 
(|Castellani et all l2006j )_ is near the average for Galac- 
tic globular clusters (jSalaris et al . 2004). However, this 
ratio could retain a "normal" value if giants are re- 
moved exclusively from the upper RGB (which is the 
most sparsely populated), or if the stars avoiding the 
HB phase (by never igniting He) are small in number. If 
there is a large population of He white dwarfs, they may 
be detectable at the bright end of the white dwarf cool- 
ing sequence because they cool more slowly than the CO 
white dwarfs produced b y traditional e v olutio nary paths 
(jCastellani et al.ll2OO60 . iDieball et all (|2005f) identified 
about 40 white dwarf candidates using STIS UV imagery, 
and found the numbers to be in rough agreemen t with 
theoretical predictions. Unfortunately, the Dicball et al.l 
survey does not reach far d own the WD sequence, and 
models (jSerenelli et al.l[2002l ) predict that the cooling age 
of a He WD is only about twice the cooling age of a CO 
WD near their completeness limit. (The ratio becomes 
much more extreme the fainter one goes.) Thus, if the 
He WD population is less than the size of the observed 
population, the production rate for He WD must be less 
than about half the CO WD production rate. 

If the true luminosity function of RGB stars is given by 
the Victoria-Regina or Teramo models, then NGC 2808 
shows a deficit in logA^ of at most about 0.1 (20% in 
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TV) in the magnitude range where He WDs are proba- 
bly produced (J — It rob ^ 0.4). If these "missing" gi- 
ants produce He WDs directly, then this would enhance 
the observable population of WDs in the iDieball et al.l 
sample by about 40%. The uncertainties in the theoreti- 
cal predictions and in the WD numbers are not yet able 
to rule this possibility out. A slightly deeper survey of 
the cluster's WD population would help definitively set- 
tle whether there is a significant populatio n of He WDs. 
Such a survey has been done for u) Cen (jMonelli et al.l 
|2005() . another cluster with a extensive blue HB tail, and 
again the results are in rou gh agreement with evolution- 
ary timescales. However, iMonelli et al.l also note that 
the observed WDs appear to be redder than expected, 
and one of the possible explanations of this is a large 
population of He core WDs. 

4. DISCUSSION 

Because our analysis implies that the RGB evolution 
is virtually independent of chemical composition and 
age inputs, it should be possible to merge photomet- 
ric data for stars from clusters with relatively hetero- 
geneous characteristics. Multiple cluster samples could 



place very tight constraints on physics inputs like neu- 
trino losses (and non-standard neutrino emission mech- 
anisms as well). Indeed, our comparisons indicate that 
the neutrino losses used by the Yale-Yonsei models can 
already be ruled out based on the comparisons with the 
massive globular clusters M5 and NGC 2808. 

Although our method can only be applied in a practical 
way to the most massive globular clusters, it does provide 
a new way of probing the "second parameter" problem in 
horizontal branch stars. If there is a dynamical influence 
on cluster giants causing them to leave the RGB early 
in NGC 2808, then we would expect similar features to 
be present in other clusters with extreme blue HB tails. 
Alternately, for a cluster like 47 Tuc with no obvious blue 
HB extension, the upper RGB should match models. 
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Fig. 1. — Comparisons of theoretical cumulative luminosity functions after shifting in magnitude to match the tip of the RGB in I 
band, and normalizing just above the RGB bump. The models are Victoria -Regina (VandeiiBerg et al. 200(j), Teramo (Cassisi ct al. 2004), 
Yale-Yonsei liPemargue et al.ii2004i ). and Revised Yale H Green et al.|[T98^ isochrones. In the lower right panel, the Victoria-Regina and 
Teramo models overlap. The two Yale-Yonsei models with T^f ^-color transformations from IGreen et al.l 119871 {short dashed line) and 
ILeieune et al.|[l998 (long dashed line) fall higher. 




Fig. 2. — The cumulative luminosity function for bright RGB stars for NGC 2808 (this work) and M5 l|Sandauist fc Boltell2004l '). along 
with a Victoria-Regina model tVandcnBcrg ct al. 2006.') for [Fe/H] = -1.14 and age 12 Gyr. 



